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AB SURACT 


MiemcecieeveiryecrucK problem iS one in which a truck is 
loaded with m packages, one package to be delivered to each 
wa Gesuinacvions. The amount of fuel consumed by the truck 
is directly dependent upon the Sieur total weight of the 
truck, which includes both the weight of the packages and 
mre amount of fuel remaining in the tank. The problem is 
to determine a sequence in which to deliver all m packages 
Were will minimize total fuel consumption. A branch and 
bound algorithm for obtaining optimal solutions to the 
delivery truck problem is presented, along with several 
sample problems with their solutions. A brief report 


On computational experience is included. 
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ieee cee OUVUe ELON 


tie CeliververtrucKeproblem 1S one in which a truck is 
loaded with m packages, one meter ane to be delivered to each 
Sue cescinations. The amount of fuel consumed by the 
meek 1S directly dependent upon the current total weight 
Sechne truck, which includes both the weight of the packages 
mee cone amount of fuel remaining in the tank. The problem 
is to determine a sequence in which to deliver all m packages 
faeeen Will minimize total fuel consumption. 

ihees ron lLem @s “closely related to eats de reyes ll ab ates 
Salesman problem. It differs in that the route is not 
eeteciea and that fuel consumed en each leg of the journey 
imommoc Only dependent upon the length of the leg, but also 
mmemcurrenl. weight of the truck. 

This thesis is the presentation of an algorithm for 
eeeeining optimal solutions to the delivery truck problem. 
The methodology of breamech vanc bound ts used to implicitly 
enumerate all possible solutions. Several example problems 
with solutions are presented followed by a brief report 
of computational experience for problems with m ranging 


from five to twenty. 





ae OD) RESEARCH 


The delivery truck problem is clearly a problem in 
miesarea Of discrete mathematical optimization. Reference 
8 is a current summary of results applicable to such 
aren lems. 

PUmeciemrescarci on tie erdinary traveling salesman 
eaeaprem is Summarized in Ref. 2, in which, after an analysis 
PemcomputraviOnal results, it is recommended that dynamic 
programming [Ref. 5] be used for problems with thirteen 
cities or less, and that branch and bound [{Ref. 9] be used 
for the rest. | 

Pimlico lene Of OOLaIning Gpvuimal solutions to the 
feeuriea traveling salesman problem in which no return to 
miemwsatarting city 1s required is presented in Ref. /. 

In an extension of the traveling salesman problem called 
"The Delivery Problem" [Ref. 4] more than one salesman is 
available to eeu Tie Dredecermanied Sel 701 CiC1eS. 

‘Another Cxtvension of the traveling salesman problem, 
sometimes referred to as the "Milk Bottle Problem", [Ref. 4] 
meeonay in which the milkman, (salesman), is to deliver 
Pier oGiOnempouvr eso. milk to a known set of houses, 
(cities). Hewimnust, Walk and wishes to minimize the amount 
of work expended carrying the milk bottles. The problem 
imc i meosei ere laved Go one Gelivery truck probien, 


Picmaitiuercemice being thal mo account 1s taken for fuel 





used along the way. A branch and bound algorithm for the 
[mee DOGG Le problem 2S presenved in Ref. 3. This algorithm 
is a direct extension of the work of Little, et al., [Ref. 
& le 

The delivery truck problem can be formulated as an 
Meese ifsavion problem. Current results in Investigation 
fimeery are contained in Ref. 1. The algorithm presented 
iaeenis thesis uses an approach similar to that of the 


Dranch and bound algorithm presented in Ref. 1. 





Tit. MATHEMATICAL DESCRIPTION 


ihestORlOWwire mMotalion will be used. Let 


Can.) =e Of GClstances between destination 


iJ 
1 omemecsuBiau On |. 31. = 1 2 yee jn, 
i Se erGe eran erae ie 
TT = (mr, ; Toaeeey T,) be a permutation of 
Miewiitee lero tc wen. Pepresentang an 
Oma mineewinien Gdesi-inat LON Sd 2 o..4 5 tl 
ByaS  NalSligs ol: 
t = the amount of fuel required al niss eye ho) omen 
foro llowed. 
1 = the amount of fuel required if destinations 
Ts To, USD Te Ts are-vVisived in -vuhat 


order... FCleariy {| = =f 
T il 


fuel required to travel from destination 


‘xj 
on 
tte 
Lama) 
i 
Sw 
i} 


Tr; UGr Gesvimation ue given that i destina- 
Tions have already been visited and fe 


al 
fuel has been consumed. 


Let fs = QO. Then 


f= & F Giese os 
=( Une gle 


MicemeOc livery muecK DprROplem 15 to select that path mt which 


minimizes ties 








tae Sur oe whan CoOnvaimed herein assumes it is known 
foeen destination will be the first visited, and the problem 
mmerps Witch the truck initially located there. If the 
mieest destination is not known, the truck can be assumed 
to be located at a dummy location added to the given set 
Sumaestinations, with travel allowed to any valid destination 
eueezeroO fuel consumption. 

Given that the truck is initially located at destination 
one, there are (n-1)! possible solutions. 

imam rerenc variations of the delivery truck problem, 
eerrerent fuel consumption functions F are appropriate. In 
Memeee thesis, F is taken to have the following form: 


: z k 
F Gi, i .) = cc, = ic. - i Wel 


ie re 
Cc) .= an Gebit~arrly wenMesen constant 
Cy = PRColLmbhiy soeodegewelphe sot the truck 
Kee. = 
oS —MGemVewoncmwoer mackag@e Ki kK = 1,...4,n. 
NOGer sr ate. t ee Smee iirc DreoLen reduces 
iL 
Peco mil~eborele Hroblem of Ref. 3. If, in addition, 
co; 0 i weieeune problem reduces to the modified 


uraveline salesman problem in which no return to the 


weet eC lL) do. PeCULPed. 











eure Se 


For the experimental testing discussed in Section IV(E) 


the constants chosen were c, = SO ea = 8000 and 
8 ie a Ore eae ec. oN ne resulvane function 7s 


MeececOo illustrate variation in fuel consumption due to 
emer loading as shown in Figure 3.1. bine a represents 
consumption Cie Orand 1 = 0, sbine b. for 1 = 5, f£ = 0, 


emoeline ¢ for i = 5 and f = 500. 





ite oe AGORA 


Mies seccriom besins by describing how the branch and 
bound method is applied to the delivery truck problem, 
followed by a general statement of the algorithm. (A flow 
eeeoeram Ol Che algorithm is contained in Appendix A and 
the computer program immediately afterwards.) Several 
faavere problems are used to illustrate the output of the 
ro rivimn. 2 Oriel summary Of computational experience 
is presented. 

A. APPLICATION OF BRANCH AND BOUND TO THE DELIVERY 
TRUCK PROBLEM 

The branch and bound method consists of sélectively 
Mercictioning the set of all feasible solutions and computing 
Beounas On the objective function for each element of the 
meee on. This process is continued until one element 
memene Parcition containing a single solution is obtained 
for which the associated bound is at least as good as that 
of any other element. 

A branching tree 1S initiated and extended such that 
a Tmeeode mate trineweaag Of a branch of the tree corresponds 
fora SEL Of SOlUtIONS. The set of all such nodes specifies 
a partition of the solution space. The bounds associated 
aimed chimetlemenymor The sDartacrion, and therefore with 


eJeremove wmmnncmciadsc! a branch of the tree, represent a 


10 





mewer bound on the amount of fuel used if any solution 
fenvermea in that element of the partition is followed. 
Associated with every node of the tree is a destination. 

Peet olal Wode, representing all soluticns, corresponds 
memcdestination number one, the initial location of the 
mance. Branching corresponds to selecting a node from 
Mimeeh: to branch and in specifying a destination for use 
Mmmexcenaing the pavh Specified by that branch. The branch 
mere sponding tO node one naturally is selected initially 
Mmemmexvcnision. This branch is extended to the right and 
Preeti vcomby adaqinge vwo new nodes to the tree, See Figure 


4.1. The node labeled k represents all solutions in 


Figure 4.1 Initial branching 


which the truck starts at destination one and proceeds first 
to destination k. The node labeled -k represents all solutions 
iinet rie mtbruck Goes not @O to destination k first. 


daa 





Associated with each node at the end of a branch of 
Mme cree 1S a path, starting with destination one, and 
me@eliding each destination number corresponding to a 
positively labeled node in the order in which they appear 
iweeome branch. 

Fiehmomoledge Omer Neus rOoetomGe  eclred forsexvension ,; 
me path specified by that branch is called the current 
path. 

The bound on each ending node with a positive label 
applies to all solutions starting with the path specified 
Dy that node. The bound on those with negative labels 
corresponds to all solutions starting with the path specified 
which do not include as their next element, all negative 
Miewews Ol nodes at the end of that branch. In illustration 


of this, consider Figure 4.2. The bound on node labeled 5 


Figure 4.2 Sample tree 


172 





meee s Fo all solutions starting with elements 1, 3, 5, in 
wees order. The bound on the node labeled -5 applies to 
PeeecoOlUCiOnS Starting With elements 1, 3, and not having 
fmectimomas the next element. 

Mower oraien 1s Seleevedq tor Extension, it is necessary 
Mmemecebermine which destinations are eligible for use in 
epevending the current path. This is done by determining 
Pemeecestinations not in the current path and not prohibited 
Mieco negatively labeled nodes at the end of the branch. 

mye oraer of an element of the partition reduces to 
Smee Wien there are no destinations eligible for use in 
eeerenading the branch corresponding to that element. When 
such an element is located with a bound less than or equal 
memene bound on all others, the path corresponding to the 


meee SOlUtION contained in that element is optimal. 


Be GENERAL STATEMENT OF THE ALGORITHM 
Bpecitiie methods of branching and bounding are ignored 
meomgerarily while the general structure of the algorithm 


is presented. 


Pabeori hint: 

Mec reave tne initial node corresponding to destination 
number one. Compute a bound for this node. Go to 2. 
aoe lCeoumaseranca Fo extend. This braneh determines the 
SrerehwmoalhsaaeGo tO 3. 

3. Determine the set BE of destinations eligible for use 


imeextending the current path. Put into E all destinations 


IES 





not in the current path. Remove from E all indices corres- 
Pending to negatively labeled nodes at the end of the branch 
eames extended. If |F| = 0, stop. If not, go to 4. 

4, Select from E a destination k for use in extending the 
Gumrent path. Extend the branch to the left and right 
Sreatang two new nodes labeled k and -k. Go to 5. 

5. Compute bounds for branches ending with nodes labeled 


k and -k. Go to e. 
oem vcCrmrnabion vne current path 1S optimal. 


C. BRANCHING AND BOUNDING RULES 
ie) Sranching 

Dranchin@etSeareWwo Part Operation, consisting first 
@eemscelecting a branch to extend, and next of selecting a 
Peosuinavion for use in extending that branch. It has 
proved to be advantageous in branch and bound to use as 
the first of these steps the selection of that branch 
whose associated bound is smallest and specifying tie 
Mmeeemeine rules in the event they occur. This rule is used 
ieee alsorithm. 

The second part of the branching operation selects 
a destination for use in extending the current path. The 
algorithm selects that eligible destination which is closest 
Menune last element in the current path, or, in other words, 
the Unvisited destination closest to the truck's current 


DOs)>FLone 


14 





e. Bounding 


The bound associated with each branch specifies a 
lower bound on the amount of fuel used if the path specified 
Memethat branch is followed. In illustration of the method 
Meea, Suppose branching is being carried out from a node 
Paieewed j CoO destination k, and that there are r destinations 
in the current path. The exact amount of fuel required to 


' follow the current path is f,, and the fuel required to 


gi 


travel from destination j to destination k is Fy (r, fi). 


The lower bound on the node labeled k is 


f, 4- Br, f 5) + L 
meere lL is a lower bound on the fuel required to visit the 
remaining n-r-l destinations. 

ei smomcolmnoo Dy PeCOonmortocrane the Upper triangle of 
the reduced matrix (d. 5) where i and j are not destinations 
mame CUrrent path, selecting the Smallest element in 
row k and the remaining n-r-2 smallest elements, ordering 
ieec G'iStances in non-decreasing order, and computing the 
amount of fuel required to traverse these distances in 
the prescribed sequence. 

Only the upper triangle need be considered since 
(d,s) is symmetric and travel between destination ilewect eis! 
CooL ner womens Carri1eGg ©UL ab most once. 

It can be seen that L is a valid lower bound by 


noting that after leaving destination k, the truck must 


IU 





meaverse an additional n-r-1 legs in its required journey. 
Glearly the shortest distance it will be required to haul 
n-r-1 packages is the smallest element in row k of the 
MmecUced Matrix. A lower bound on the next leg is clearly 
the fuel required to haul n-r-2 packages the smallest 
Setance In vhe reduced matrix, excluding the element used 
Secs irsy Jee, This procedure 1S continued until all 
meee legs of the journey are accounted for. 

When the bound is computed for the node labeled -k, 
mpemprocedure 1S the same except that selection of a distance 
for the first additional leg after leaving destination j 


is not restricted to come from row k of the reduced matrix. 


eee ROOF OF OPTIMALITY 


Pao eocemune SCUlOlvalincOlUutions to the problem, 1.¢e., 


pits a2 permucacion of anvegers 1,2,...,n} , 


poromet F be a partition of S such that 


Let De be a lower bound on all solutions contained in 


Po; Cea ¢ 


DS ee POrral. 7 in Poe 
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Miesceo cc Or the algorithm a branch of the tree is 
selected for extension. Selecting a branch for extension 
@erresponds to selecting an element of the current partition 
for refinement. If element P, is selected and [P| = 1, 
Mempurcvhner refinement is possible, and in Step 3, || aan Onn 
Pee ois point all n destinations are contained in the 
weerent path, and bound Dy isemo longer a lower bound, but 
Peect. Note that the method of selection used for oe as 
described in Section IV(C)1, is to select partition element 


meeouch that 


d, = min DL ene yesiiudl Pi aie ge nae 
Hence , P haoceamime SAC mOOCUNnG Wworcm 3s Less than or equal 
womune lower bounds associated with all elements of P. 
Mae single solution tw in . Ho eNews SOlMedOr whach 


minimizes fe 


E. oAMPLE PROBLEMS AND SOLUTIONS 
ine tollowinge are sample problems with optimal solutions. 
The x and « G@GREMACER of the destinations were selected 
meom a uniform distribution over the interval (0, 10). 
Coordinates were adjusted slightly when necessary to 
insure that each pair of destinations was at least a unit 
Gesoance aparece, The truck is assumed to ‘be initially 
Pecated at destination one. The path shown is the one which 
Tio ueswcactwel the remaining destinations and uses the 


Mmaidimum amount of fuel. 
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meee COMPUTATIONAL EXPERIENCE 

Bere eavtenaleeypertvence with the algorithm is limited, 
Dut an investigation was made into the way the time required 
to solve ann destination problem varies with n. Seventeen 
sample problems were developed as described in the preceeding 
Beeston; five with n equal to five, six with n equal to 
imomeemend Six with n equal to fifteen. The average times 
required to solve each set of problems, along with the 


average number of bounds computed, is shown in Table 4.1. 


average average # 
time (sec.) of bounds 





eae Hs 1] 


Although these results are sparse, this data suggests 
Meaue came required increases exponentially with n, as seen 
in Figure 4.9. The circles show the data points of Table 
4,1. The dots show times required to solve sample problems 


meving nm equal to 20. 
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APPENDIX A 


FLOW DIAGRAM OF COMPUTER PROGRAM 


pom toLlowilNe iOlaltiom, Nee previously defined, is used: 
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ie whos Mecemere ated - 

Ceaeinartonemumoer associated with node I, 
inode mee a 4a 5 go LN 

Viewed Mincamavely trode. T dn the bDranchings 
pice Came mOron ee 22 0 stew LN 

node selected ie re Secs 1 oO! 
Cesvinaudemeseleeved for use in extending 
the current path 


bound on node I 


Prem = Beit) 1s mode [Ets at the end of a branch in the 


NPL 


MN 


Pee Y ,Z ) 


Bree, and equals M olnerwase, where M is 
Ienateters 

MuUMO Cie wCteslaemares Cumrent ly delivered 
destination SaewoLeh the truck is currently 
located 

fuel required to travel distance Z if X 
packages have currently been delivered and 


Y fuel consumed doing so. 
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